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Abstract 

A Matlab program is presented that computes derivative corrections in the S- 
dual invariant formulation for IIB graviton scattering to any order in perturbation 
theory. The coefficients of the four-point function are produced, pertaining to the 
non-logarithmic terms. The program can be modified to find coefficients of the higher- 
point functions. Instantons have not explicitly been included. 



The graviton scattering amplitude in IIB superstring theory has been investigated 
recently in the context of manifesting S-duality. This manifestly S-dual graviton 
scattering has been analyzed systematically in pQ, and references therein. 

This text contains a computing program that automates the perturbative calcu- 
lations. To any order in derivatives this program computes the analytic portion of the 
graviton scattering, in the flat ten-dimensional superstring theory. The D-instantons 
are allowed to be included, but their form is not placed in the program. 

The program is written in Matlab, and follows the procedure in pQ by partitioning 
numbers and attaching weighted trees to these numbers. 

Given the large number of papers on the topic, the author decided to place it on 
the archive. The Matlab program may be found by downloading the TeX source. 

The quantum gravitational action may be numerically studied with this program, 
at the D fc i? 4 level, with any set of integers k. 
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Matlab Program 
clear all 

Nmin=l; Nmax=10; 

trialnumber=200; 

yp=0; flagno=l; 

nullcase=0; 

trialnumber=200 

for i=Nmin:Nmax 

for j=l:round(Nmax/3) 

Oone(i)=0; Partitiontype(i,j)=0; 

end 

end 

for N=l:Nmax 

for Np=l:round(Nmax/3); 

x(l,l)=0; p(l,l)=0; rp(l,l)=0; flagnumber(l)=0; 

nullcase=0; a=0; dimp=0; flagno=0; timem=0; 

totalfiagno=0; v(l)=0; yp=0; 

clear x,p,rp,flagnumber,nullcase; 

clear a,dimp,flagno; 

clear timem,totalflagno,v,yp; 

clear i 

if Npj=N 

x(l,l)=0; clear x p(l,l)=0; clear p 

for j=l:trialnumber 

for i=l:Np 

x(ij)=0; p(ij)=0; 

end 

end 
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for j=l:trialnumber 

flagno=l; timem=0; nullcase=0; 

while flagno==l 

yp=0; 

for i=l:Np 

x(i,j)=2*round(((N-l)/2-l)*rand)+3; 

yp=yp+x(i,j); 

end 

flagno=l; 
if yp==N 

P(:j)=x(:j)l 
flagno=0; 

end 

timem=timem+ 1 ; 

if timem^trialnumber 

flagno=0; 

nullcase=l; 

end 

if nullcase==l 
sp=0; 

j=101; 
end 
end 
end 

if nullcase =1 
for j=l:trialnumber 
for m=l:trialnumber 
for a=l:Np 
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for b=l:Np 
if ajb 

if P(aj)ip(bj) 

testl=p(a,j); test2=p(bj); 

p(bj)=testl; p(a,j)=test2; 

end 

end 

end 

end 

end 

end 

v=size(p); dimp=v(2); 
for i=l:dimp 
flagnumb er ( i ) = ; 
end 

for i=l:dimp 
for m=l:dimp 
if m==i 

flagnumber (m) = 1 ; 
end 

if m =i 

if p(:,m)==p(:,i) 
if iim 

flagnumb er ( m) = ; 

else flagnumber(m)=l; 

end 

end 

end 
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end 
end 

totalflagno=0; 
for i=l:dimp 

tot alflagno = t ot alflagno +flagnumb er ( i ) ; 
end 

rp(l,l)=0; clear rp 

for i=l:totalflagno 

for j=l:Np 

rp(j,i)=0; 

end 

end 

sp =0; k=0; 

for i=l:dimp 

if flagnumber(i)==l 

k=k+l; 

rp(:,k)=p(:,i); 

sp=sp+l; 

end 

end 

sp2=sp; test=0; 

for k=l:sp2 

for 1=1 :Np 

test=test+rp(l,k); 

end 

end 

sp=test/N 
end 
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end 

for i=l:Np 

Partitiontype(l,i)=0; Partitiontype(2,i)=0; 
end 

Partitiontype(N,Np)=sp; 
end 

for j=l:Np 

Oone (N) =Oone(N) +Partitiontype (N,j ) ; 
end 

clear sp; 
end 

Nmin=l; Nmax=10; 

clear i 

tau=i*l; 

taul=real(tau); tau2=imag(tau); 

trialnumber=200; 

A1=1;A2=1; 

for i=l:Nmax 

tensorc(i)=l; 

end 

for i=l:Nmax 

for j=l:round(Nmax/3) 

Amp(i,j)=0; 

end 

end 

for N=l:Nmax 

for Np=l:round(Nmax/3); 

if Partitiontype(N,Np) =0 
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yp=0; flagno=l; 
nullcase=0; 

x(l,l)=0; p(l,l)=0; rp(l,l)=0; flagnumber(l)=0; 

nullcase=0; a=0; dimp=0; flagno=0; timem=0; 

totalflagno=0; v(l)=0; yp=0; 

clear x,p,rp,flagnumber,nullcase; 

clear a,dimp,flagno; 

clear timem,totalflagno,v,yp; 

if Npj=N 

x(l,l)=0; clear x p(l,l)=0; clear p 

for j=l:trialnumber 

for i=l:Np 

x(i,j)=0; p(i,j)=0; 

end 

end 

for j=l:trialnumber 

flagno=l; timem=0; nullcase=0; 

while flagno==l 

yp=0; 

for i=l:Np 

x(i,j)=2*round(((N-l)/2-l)*rand)+3; 

yp=yp+x(i,j); 

end 

flagno=l; 
ifyp==N 

P(:j)=x(:j); 
flagno=0; 

end 
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timem=timem+ 1 ; 

if timem^trialnumber 

flagno=0; 

nullcase=l; 

end 

if nullcase==l 
sp=0; 

j=101; 
end 
end 
end 

if nullcase =1 

for j=l:trialnumber 

for m=l:trialnumber 

for a=l:Np 

for b=l:Np 

if ajb 

if p(aj)ip(bj) 

testl=p(a,j); test2=p(bj); 

p(b,j)=testl; p(a,j)=test2; 

end 

end 

end 

end 

end 

end 

v=size(p); dimp=v(2); 
for i=l:dimp 



flagnumb er (i ) = ; 
end 

for i=l:dimp 
for m=l:dimp 
if m==i 

flagnumber(m)=l; 
end 

if m =i 

if p(:,m)==p(:,i) 
if i^m 

flagnumber(m)=0; 

else flagnumber(m)=l; 

end 

end 

end 

end 

end 

totalflagno=0; 
for i=l:dimp 

tot alflagno = t ot alflagno +flagnumb er ( i ) ; 
end 

rp(l,l)=0; clear rp 
for i=l: tot alflagno 
for j=l:Np 

rp(j4)=0; 
end 

end 

sp =0; k=0; 
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for i=l:dimp 

if flagnumber(i)==l 

k=k+l; 

rp(:,k)=p(:,i); 

sp=sp+l; 

end 

end 

instanton=0; 
w=size(rp) 
dimbase2N=w(l); 
Q=0; 

AmpNNp=0; 
for m=l:w(2) 
R(m)=0; T(m)=l; 
end 

for i=l:dimbase2N 
S(i)=0; zs(i)=0; 
end 

for m=l:w(2) 

if rp(:,m) =0 

for i=l:dimbase2N 

zs(i)=rp(i,m); 

end 

for A=l:base2N 
Q=A-1; 

for l=l:dimbase2N 
if 1 =dimbase2N 

Q=P; 



end 

if l==dimbase2N 
if Q==l 
a(l)=l; 
end 

if Q==0 

a(l)=0; 

end 

end 

end 

k=N-3; 

s=3/2+(N-3)/2; 
if mod(k,2)==0 
gmax=l/2*(k+2); 
end 

if mod(k,2)==l 

gmax=l/2*(k+l); 

end 

for i=l:dimbase2N 

if a (i)==0 

if instanton==l 

S(i)=0; 

end 

end 

if a (i)==l 

if instanton==l 

S(i)=0; 

end 



end 

T(m)=T(m)*S(i); 

end 

end 

R(m)=R(m)+T(m); 

end 

end 

end 

end 

for c=l:w(2) 

Amp(N,Np)=Amp(N,Np)+R(c); 

AmpNNp=AmpNNp+R(c) ; 

end 

end 

end 

end 

for i=l:Nmax 
Coeff(i)=0; 

for j=l:round(Nmax/3) 

Coeff(i)=Coeff(i)+tensorc(i)*Amp(i,j); 

end 

end 



13 



